Supersymmetric Quantization of 
Anisotropic Scalar- Tensor Cosmologies 



James E. Lidsey* 

Astronomy Unit, School of Mathematical Sciences, 

Queen Mary and Westfield, 
Mile End Road, London El 4NS, United Kingdom 

P. Vargas Moniz ti§ 
Grupo de AsTrofisica e Cosmologia (GATC) 
Departamento de Fisica, Universidade da Beira Interior (UBI) 
Rua Marques d'Avila e Bolama, 6200 Covilha, Portugal 

February 7, 2008 



Abstract 

In this paper we show that the spatially homogeneous Bianchi type I and 
Kantowski-Sachs cosmologies derived from the Brans-Dicke theory of gravity 
admit a supersymmetric extension at the quantum level. Global symmetries in 
the effective one-dimensional actions characterize both classical and quantum 
solutions. A wide family of exact wavefunctions satisfying the supersymmet- 
ric constraints are found. A connection with quantum wormholes is briefly 
discussed. 



PACS number(s): 98.80.Hw, 04.60.Kz, 11.30.Pb 



*E-mail: jel@maths.qmw.ac.uk 
TE-mail: pmoniz@mercury.ubi.pt 



■^URL: http:/ /www.dfis.ubi.pt/^pmoniz 



Also at CENTRA, 1ST, Rua Rovisco Pais, 1049 Lisboa Codex, Portugal 



1 Introduction 



Quantum cosmology applies the fundamental principles of quantum physics to the 
entire universe. (For a review, see, e.g., Ref. [|I|]). The wavefunction of the universe is a 
functional on the configuration space (superspace) and obeys an infinite-dimensional 
partial differential equation - the Wheeler-DeWitt equation 0. In view of the severe 
technical difficulties that arise in solving this equation, the normal procedure is to 
arbitrarily confine the fields to the neighbourhood of spatial homogeneity. Effectively, 
the infinite number of inhomogeneous modes and their interactions are truncated out 
and the configuration space (minisuperspace) is therefore finite-dimensional. The 
Wheeler-DeWitt equation then determines the evolution of the wavefunction on the 
minisuperspace and a given trajectory mapped out by the wavefunction may be in- 
terpreted as a cosmological space-time. 

The validity of the minisuperspace approximation remains an open question to 
date. It is clearly inconsistent with the uncertainty principle since the amplitudes and 
momenta of the inhomogeneous modes are assumed to vanish simultaneously. Kuchar 
and Ryan addressed this question quantitatively within the context the Bianchi type 
IX cosmology and found that imposing additional symmetry on the model altered 
the nature of the physical predictions || (see also Refs. 0, On the other hand, 
Sinha and Hu employed the techniques of coarse-graining and interacting field theory 
to derive a condition that must be satisfied for the approximation to be justified ||. 
Difficulties in coupling non-trivial spinor fields to highly symmetric spacetimes have 
also been highlighted by Henneaux 0. 

Despite these uncertainties, however, the expectation is that the main features of 
the wavefunction should be preserved in a more general analysis || . In principle, the 
wavefunction of the universe yields the probability that a spatial hypersurface evolves 
from a given initial state. However, ambiguities arise when attempting to invoke such 
an interpretation due to the hyperbolic nature of the Wheeler-DeWitt equation: a 
conserved current with a positive-definite probability density is not possible. 

One possible resolution of this and related difficulties is to extend the standard 
quantization of the universe in a super symmetric fashion. (For a review, see, e.g., 
Refs. H |T0j] ) . Supersymmetry may help in the quantization of gravity for a number 



of reasons. Indeed, earlier work on supergravity theories |TT|] and recent develop- 
ments in superstring theory |TJ| and M-theory |13] indicate that supersymmetry is 
a fundamental ingredient of any unified description of the fundamental interactions. 
Consequently, an analysis of the very early universe that includes supersymmetry 
is well motivated and a number of authors have developed models of the early uni- 
verse where both quantum gravitational and supersymmetric effects are important 



H [K| [L4|, 0, |T(| [17], |18| . The advantage of such an approach is that the quantum 
state of the universe, is annihilated by supersymmetric constraints that are linear, 
first-order differential equations in the bosonic momenta variables. This is in contrast 
to nonsupersymmetric quantum cosmology, where a second-order Wheeler-DeWitt 
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equation has to be solved subject to suitable boundary conditions [O, |20], |22|]. The 



super symmetric algebra necessarily implies that \& also obeys the Hamiltonian con- 



straint and it is therefore sufficient to solve the first-order constraint equations [23 
In many cases, this resolves ambiguities in the choice of factor ordering. Further- 
more, supersymmetric quantum cosmology places the results of standard quantum 
cosmology in a wider perspective |] and a study of quantum minisuperspaces with 
supersymmetry may also provide (in spite of the obvious truncations) some helpful 
insights concerning the set of states that represent a complete formulation of quantum 
gravity with the other interactions. 

In recent years, two attractive (and possibly related) approaches to supersymmet- 
ric quantum cosmology have been developed. One approach is to begin with N — 1 
super gravity [|11], |T^] in four dimensions and reduce the system to a one-dimensional 
model by invoking a suitable homogeneous ansatz || |10|, [15], [16| . This leads to a 
minisuperspace with N = 4 local supersymmetry. Alternatively, one may integrate 
a purely bosonic action over the spatial variables to derive a (1 + 0)-dimensional la- 
grangian and then perform a supersymmetric extension of the corresponding Hamil- 
tonian system by employing the quantization rules of the supersymmetric a-model 
24, 25, 26, [17], 18]. This results in an N = 2 supersymmetry. In particular, this 
process could be related to the fact that any one-dimensional system is supersym- 
metric provided its ground state is normalizable f27| . Moreover, this technique can 



be generalized to higher dimensions by employing Darboux transformations 

In this paper we employ the latter approach to quantize spatially homogeneous 

0. The 



cosmologies PDj within the context of the Brans-Dicke theory of gravity 
Brans-Dicke theory is relevant to the early universe and arises as the effective action 
of higher-dimensional gravity theories and, in particular, superstring theory [12 . 
Moreover, the spatially flat and isotropic Brans-Dicke cosmology exhibits a discrete 
'scale factor duality' |TB|, |5T| . This symmetry forms the basis of the pre-big bang 
inflationary scenario []3l], [32|] (for a review see, e.g., Ref. and its origin can 

be traced to the T-duality of string theory |34j]. The consequences of scale factor 
duality for string quantum cosmology have been explored by a number of authors 
|18|, [35|. In particular, supersymmetric quantum states have been found that respect 



the duality symmetry of the classical Hamiltonian |TJ|. Our purpose in this paper is 
to perform a supersymmetric extension of more general spatially homogeneous cos- 
mologies. Specifically, we consider the locally rotationally symmetric (LRS) Bianchi 
type I model and the Kantowski-Sachs universe [29]. Generalizations of scale factor 
duality have been shown to exist in these models [35, |57|] and we find supersymmetric 
wavefunctions that respect these symmetries. 

The paper is organized as follows. In Section II, the global symmetries of the ac- 
tions and the supersymmetric quantization procedure are reviewed. The LRS Bianchi 
type I model is quantized in Section III and the vacuum Kantowski-Sachs model is 
quantized in Section IV. We conclude with a discussion in Section V, where we also 
comment on the possible relationship between the supersymmetric minisuperspace ex- 
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tension [17, [18|, ^7], 28| and minisuperspaces retrieved from more general supergravity 
theories g, JTOj, [OJ. 

We assume throughout that h = 1. 



2 Supersymmetric Quantum Bianchi Cosmology 



2.1 Duality and the Wheeler— DeWitt Equation 

We consider the four-dimensional Brans-Dicke action given by 



S 



d Xy 



R - w(V$) 2 - 2A 



(2.1; 



where R is the Ricci curvature of the spacetime with metric g^ u , g = detg^ u and 
$ represents the Brans-Dicke (dilaton) field. The coupling between the scalar and 
tensor fields is parametrized by the constant, u, and A is the cosmological constant in 
the gravitational sector of the theory. A consistent truncation of the string effective 
action is given by Eq. fl2.ip for uo = — 1 and A < [12|. Dimensional reduction of 
higher-dimensional Einstein gravity on an isotropic, ci-dimensional torus results in 
the above action, where uo = — 1 + 1/d and $ determines the volume of the internal 
space rC> 
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The metric for the class of spatially homogeneous, LRS cosmological models with 
constant time hypersurfaces containing two-dimensional surfaces of constant curva- 
ture, k, is given by 



ds 2 



-N 2 dt 2 + e 2a ~^dr 2 + e 2a+2/3 d£l 
-N 2 dt 2 + a\dr 2 + a\d£l\ k , 



2 

2,k 



(2.2) 
(2.3) 



where N is the non-dynamical lapse function, dQ% k is the unit metric on the constant 
curvature two-surfaces, e 30 ^ = aia 2 determines the effective spatial volume of the 
universe and (3 = (l/3)[ln(a 2 /ai)] determines the anisotropy of the model. The cases 
k = { — 1, 0, +1} correspond to the Bianchi type III, I and Kantowski-Sachs universes, 
respectively. The geometry of the spatial sections of the Kantowsk-Sachs model is 
S 1 x S 2 . The symmetry group of these surfaces is of the Bianchi type IX, but only 
acts transitively on two-dimensional surfaces that foliate the three-space. 

Integrating over the spatial variables in Eq. (|2.1| ) for the metric ansatz (|2.2j ) yields 
the minisuperspace action: 



S 



dtNe 



a 2 d<i> <I> 2 3 2 , 0n 9fl 



Introducing the new variables 1 



a 



'3 + 2cj 
~4 + 3u 



($ - 3a) 



(2.4) 



(2.5) 



1 We assume throughout this paper that u> > —4/3. 
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II 



'8 + 6a; 
2 + a; V ! 
1 



L-[a + (l+u;)$]+/^ (2.6) 
a + (l+u;)$-2/3] (2.7) 



v / 2T^ 

implies that we may diagonalise the kinetic sector of the reduced action (|2 



+ 2Nke {c - K)a ~ Gu - 2NAe~ Ka } , (2i 



where 



/4 + 3a; , , 

« = t/:TT^7 ( 2 -9) 



3 + 2a; 
(3 + 2a;)(4 + 3a;) 



^ ■ 2(1 + " } (2.10) 



CeJ^. (2.11) 
V 4 + 3a; V ; 

Global symmetries in these models, corresponding to a generalization of scale 
factor duality were uncovered in Ref. |36|. The action (j2.8|) is invariant under the 
discrete Z 2 'duality' symmetry 

u = u, v = —v, a = a (2-12) 

and in terms of the original variables in Eq. (|2.4j), this is equivalent to 



4 + 3a; 2(1 + a;) 4 . 

a = — ra- + t0 2.13 

3(2 + a;) 3(2 + a;) 3(2 + K 1 

$ = — a — $ + — (2.14) 

2 + a; 2 + a; 2 + uj v ; 

5 2 2(1 + a/) 2 + 3a; _ , n 1 c , 

/3 = — r« + + ^ + —, t0- (2-15) 

H 3 2 + a; 3 2 + a; 3(2 + u) H V ' 



The scale factors transform such that 



oj = a^e -- (2.16) 
a 2 = a 2 . (2-17) 

Thus, the scale factor a 2 is invariant under the symmetry transformation, whereas ai 
undergoes a direct inversion for the string inspired case, u = — 1. 
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The spatially flat Bianchi type I model (k = 0) also exhibits a global SO (2) 
symmetry that acts non-trivially on the variables {u,v}: 



u = cos 9u — sin 9v 
v = sm9u + cos9v, 



(2.18) 



where 9 is a constant. The equivalent transformations on the scale factors and dilaton 
field were presented in Ref. [|36[| . The variable, a, transforms as a singlet under Eq. 

The field equations for these models can be expressed in the form of an uncon- 
strained Hamiltonian system, where the Hamiltonian vanishes. The momenta conju- 
gate to the variables {u, v , a} are given by 



7T t . 
7T„ 



2iie- Ka 
2ve~ Ka 
-2de~ K \ 



from which the classical Hamiltonian constraint follows: 

H = -ti 2 u - n 2 + nl + 8ke {c - 2K)a e- Gu - 8Ae- 2K(7 . 
Eq. ( |2.21| ) may be written in the more compact form 

H = G ab n a 7t b + W(q a ), 



W(q a 



-Gu 



+ 8Ae 



— 2kct 



(2.19) 
(2.20) 



(2.21) 



(2.22) 
(2.23) 



where q a = (a,u,v) (a = 0,1,2) and G ab = diag(— 1, 1, 1) is the minisuperspace 
metric. By identifying the conjugate momenta with the operators ir q a = 7i a = —id/dq a 
and neglecting ambiguities that arise in the factor ordering, we arrive at the Wheeler- 
DeWitt equation: 



da 2 du 2 dv 2 



2.2 Supersymmetric Quantum Cosmology 



^ = 0. 



(2.24) 



In this subsection we summarize the procedure for attaining a supersymmetric 
extension of the models. In general, such an extension of the system is possible 
if a solution, / = I(q a ), can be found to the Euclidean Hamilton- Jacobi equation 



24, 25, 17 



G 



ab 



dl dl 



W{q a ) 



dq a dq b 

A quantum Hamiltonian, H, may be defined by the conditions 



(2.25) 



2H=[Q,Q] + , Q 2 = Q 2 







(2.26) 
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and 

[F,Q]_ = [iJ,Q]_ = 0, (2.27) 

where Q is a non-Hermitian supercharge and Q is its adjoint. The functional forms 
of these supercharges are 

Q = r + i^j (2.28) 

and 

Q = r (n a - i^j , (2.29) 

respectively, where the corresponding fermionic (Grassmannian) variables are defined 
by 

$ a = e a , V 6 = G ab 4r (2.30) 
r r d6 a y ' 

^> fe + ^V a = 0, ^> 6 + ^ a = 0. (2.31) 

Eqs. (|2.26| ) and (|2.27| ) represent the algebra for a N = 2 supersymmetry. For 
the three-dimensional minisuperspace that we are considering, the supersymmetric 
wavefunction can be expanded in terms of the Grassmann variables 6 a : 

= A + + B a 9 a + ^e abc C c 9 a 9 b + yL.fl'W, (2.32) 

where the bosonic variables {A + , B a ,C c , A_} are functions of the minisuperspace 
variables, e a & c is totally antisymmetric on all its indices and e i2 = +1, etc. The 
supersymmetric wavefunction is then annihilated by the supercharges: 

= (2.33) 
= (2.34) 

and automatically satisfies the Hamiltonian constraint due to Eq. ( p. 26 ). 

The Euclidean Hamilton-Jacobi equation for the LRS Bianchi models we are 
considering is given by 

" (§)' + (I)' + © 2 = - 8te< °- 2 ^-°" + 8A ^' < 2 - 35 > 

Thus, the problem of quantizing these models in a supersymmetric fashion involves 
finding a solution to Eq. ( |2.35| ) and then solving the simultaneous constraints Q2.33D 
and ( [2.34|) subject to the ansatz Q2.32j ). Furthermore, since the models exhibit a 
classical, global symmetry, it is natural to consider those solutions to Eq. ( |2.25| ) that 
respect this symmetry. 

In the following Section we employ this method to quantize the LRS Bianchi I 
cosmology. 
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3 Supersymmetric LRS Bianchi I Quantum Cos- 
mology 

The Euclidean Hamilton-Jacobi equation for the LRS Bianchi type I cosmology 



is 



da 



+ 



'd£ 
du 



+ 



'or 

dv 



8Ae 



-2ku 



(3.i; 



A solution to Eq. fl3~T| ) that respects the global symmetry ( |2.18|) and discrete Z 2 
symmetry (|2.12| ) of the reduced action (|2.8|) is given by 



A < 
A > 



I = T 
I = T 



VA 2 - 8Ax 2 - Acotanh 1 



vM 2 - 8Ax 2 - Atanh" 1 



'y/A 2 -8Ax r 
'VA 2 -8Ax 2 ' 



+ AVu 2 + v{3.2) 



+ AVu 2 + v 2 , (3.3) 



where A is an arbitrary constant and x = e Kcr . In the A > case one also requires 
a < ln(2v / 2A 1 / 2 /A) if A > and a > ln(-2 v / 2A 1 / 2 /A) for A < 0. 

Given the solution (|3T2"|), we could in principle quantise the system in a manifestly 
supersymmetric fashion. For simplicity, however, we consider the case 2 where A = 0. 
The Euclidean action (pj.2|) then simplifies to 



K 



and it follows immediately from Eq. (|3.4j) that we require A 
The supercharges ( [2.28D and ( |2.29| ) are then given by 



(3.4) 

-A < for consistency. 



and 



. d d . d d . d d r— _ na d 



da ou ov 



(3.5) 



(3.6) 



respectively. 

The constraint (|2.34j ) yields the set of coupled, first-order partial differential equa- 
tions 



.OA, 



— 1- 



da 



=F i2\/2\e- Ka A A 



0. 



(3.7) 



2 This corresponds to the limit a —>■ —00 and denotes a weak coupling regime for a. In the 
strong limit (a — > +00), the term in A dominates and A can be positive. In the former situation, 
the 'averaged' scale factor volume (represented by a) and the dilaton are more important, while in 
the latter a large anisotropy in the spatial directions dominates. This allows A to be positive. In 
particular, one may have a < hoo and u — » —00, with, e.g., the singularity 02 — > 0, a\ — > 00. 
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— I- 



,dA+ 

du 
.dA + 

dv 



0. 
0. 



.dB 1 
l da 


+ 


.dB 
du 


.dB 2 




.dB _ 


'da 


+ 


dv 


.dB 2 


+ 


.dB 1 


du 


dv 


.ldC° 




.ldC 


% 2 da 




% 2 du 



=F i2\/2\e~ K(T B 1 = 0, 
=F i2V2Xe- Ka B 2 = 0, 



0. 



ldC< 



— i— 



2 dv 



=F iV2\e- Ka C° = 0. 



From the corresponding constraint ( ^.33 ), it follows that 
.dA- 



' da 
.dA + 

du 
dA_ 



—% 



dv 
idC 1 

~2~da 

dC 2 



=F z2^2Ae _KCT A 

= 0, 

= 0, 
dC° 



0. 



± i V2\e~ Ka C 1 



0. 



da 
.dC 2 

du 
dB° 



du 

i-rr— =F i2V2\e~ Ka C 2 = 0, 



dv 
dC 1 



0. 



dv 

.dB 1 3B 2 



— %- 



T i2^2\e- Ka B° 



da du dv 

Eqs. (|3.7| ), ( p.8|) and (|3~^ ) immediately imply that 

A + = A° +e f 

with A° + is an arbitrary constant and 

f = 



0. 



2V2A 



Similarly, we deduce from Eqs. ( |3.14| ), ( |3.15| ) and (|3.16|) that 

A_ = A°_e~f, 
where A°_ is a second, arbitrary constant. 



(3.8) 
(3.9) 
(3.10) 
(3.11) 
(3.12) 
(3.13) 

(3.14) 
(3.15) 
(3.16) 
(3.17) 
(3.18) 
(3.19) 
(3.20) 

(3.21) 

(3.22) 
(3.23) 
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To proceed in solving Eqs. ( |3.10D , (|3.11|) , ( p,12|) and ( |3.20|) , it proves convenient 
to redefine the functions B a as follows: 



B„ 



B a e f , a =(0,1,2). 



(3.24) 



From the definition of f{cr) given in Eq. fl3.22|) , it follows from Eqs. ( [3.10|) , (|3.11| ), 
( gTTj ), ( pOD and (ggg ) that 



dB x 


dB 


da 


du 


dB 2 


dB 


da 


dv 


dB 2 


dBx 


du 


dv 


dB 


d^ 


da 


du 



0. 

o, 

dB 2 

dv 



2KfB = 0. 



Similarly, by introducing the new set of variables C b defined by 

C b = C b e- f , 



(3.25) 
(3.26) 
(3.27) 
(3.28) 

(3.29) 



we derive a new set of equations that are equivalent to Eqs. ( |3.13| ), Q3.17D , (|3.18| ) and 



dC 1 


+ 


dC° 


da 


du 


dC 2 


+ 


dC° 


da 


dv 


dC 2 




dC 1 


du 




dv 


dC° 




dC 1 


da 


+ 


du 



o. 

o, 
o, 

d& 

dv 



+ 2k/C = 0. 



(3.30) 
(3.31) 
(3.32) 
(3.33) 



By manipulating Eqs. ( |3.25| )-( p.28j ) we arrive at the following set of decoupled 
equations^] 



d 2 B 2 
da 2 

d 2 Bp 
da 2 

d 2 B x 
da 2 



2«/ 
2«/ 
2k/ 



dB 2 
da 

dBp 
da 

dBi 

~da~ 



d 2 B 2 
du 2 

2k 2 JBo 



d 2 B 2 _ 
dv 2 
d 2 Bo 



du 2 
d 2 B x d 2 B x 



d 2 Bp 
dv 2 



du 2 



dv 2 



0. 



(3.34) 
(3.35) 
(3.36) 



3 For example, Eq . ([3.34 ) is derived by applying the differential operator d/dv on E q. (|3.28| ), 
then acting on Eq. (|3.26 ) with d/da and on Eq. (3.27) with d/du. By employing Eq. ( 3.26| ), we 
then arrive at Eq. (3.34) above. A similar procedure leads to Eqs. ( |3.35| ) and ( |3.36| ). 
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Applying an equivalent technique to Eqs. ( |3.30|) -( |3.33|) results in a set of decoupled 
equations for the amplitudes C c \ 



Q2 C 2 dC 2 Q2 C 2 Q2 C 2 

2K f^r~ + + = °> ( 3 - 37 ) 



da 2 da dv 2 dv 

2K/ «? + ^ + ^ = , (3 . 39) 



da 2 da du 2 dv 

Eqs. ( |3.25|) - (|3.28|) can be solved analytically if B 12 are independent of a. Eqs. 
( |3.34| ) and (|3.36|) then imply that these variables satisfy the two-dimensional Laplace 
equation, subject to the integrability condition ( |3.27|) . Eqs. ( |3.25| ) and ( |3.26| ) further 



imply that B is independent of {u, v} and consistency between Eqs. ( |3.28|) and ( |3.35|) 
results in a further integrability constraint 

The functional form of B follows immediately up on integration of Eq. Q3.28] ) , B = 
er^ ' . It is interesting that this is also the wavefunction ( |3.23|) for the filled fermion 



sector. Similar conclusions follow for the functions C c . If C 1,2 are independent of a, 
satisfy the two-dimensional Laplace equation, the integrability condition, dC 1 /du = 
—dC 2 /dv, and Eq. (|3.32|) , then the function C° is given by the wavefunction ( [3.21D 
for the empty fermion sector, C° = eA 

Finally, it is interesting to compare the wavefunction ( |3.21| ) for the empty fermion 



sector with the general solution to the bosonic Wheeler-DeWitt equation ( 2.24 ). 
When the wavefunction depends only on the variable a, the general solution to Eq. 
( f2.24| ) is given by 

* = c 1 I {f) + c 2 K (f), (3.41) 

where 1$ and K are modified Bessel functions of the first and second kind with order 
zero, / is defined in Eq. ( |3.22| ) and q are arbitrary constants. In the large argument 
limit, the modified Bessel function of the first kind asymptotes to the form I oc 
/ _1 / 2 exp(/) and, consequently, there is a correlation, up to a negligible prefactor, with 



the fully bosonic component ( [3.21 ) of the supersymmetric wavefunction. Indeed, the 



solution A + = exp(/) is an exact solution to the bosonic Wheeler-DeWitt equation if 
a suitable choice of factor ordering is made when identifying the momentum operator 
conjugate to the variable a. In general, the ambiguity in the factor ordering can be 
accounted for |19[ by identifying 



K = - e ~ p V eP V ( 3 - 42 ) 

da da 

for some constant p in the classical Hamiltonian ( [2.21| ). In this case, the corresponding 
Wheeler-DeWitt equation is then solved by Eq. (|3.21|) for p = k. 
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4 Supersymmetric Kant owski— Sachs Quantum Cos- 
mology 

In this Section, we consider the supersymmetric quantization of the vacuum 
Kantowski-Sachs, Brans-Dicke cosmology where A = 0. The Wheeler-DeWitt and 
Euclidean Hamilton- Jacobi equations are given by 



d 2 



+ 



d 2 



+ 



d 2 



da 2 du 2 dv 



and 



'dl_ 
da 



du 



+ 8e J 



'd£ 
dv 



-Bu 



8e 



V = 



Aa+Bh 



(4.43) 



(4.44) 



respectively, where A = C — 2k and B = —G. 

We now assume?] the wavefunction does not depend on the variable v and introduce 
'null' variables over the reduced (1 + l)-dimensional minisuperspace: 



S 



s = 



A 2 



^ eXp 



exp 



-{A + B){a + u) 



-(A-B)(a-u) 



(4.45) 



The Wheeler-DeWitt equation ( |4.43| ) transforms into the unit-mass Klein-Gordon 
equation 

" d 2 



dsdr 



- 1 



* = 



and particular solutions to Eq. (|4.46|) are given by 



-ifis+ir I fl 



(4.46) 



(4.47) 



where [i is an arbitrary, complex constant. If Im/z < 0, the modulus of the wavefunc- 
tion is square-integrable. The general solution to Eq. (|4.46| ) may be expanded as a 
linear superposition of the family of solutions ( f4.47| ): 



gen 



d 2 (xF(fx,fx*)^^ 



(4.48) 



The function, F, represents a weighting function. If this is finite and only supported 
over a closed area of the Im/z < sector of the complex /x-plane, Cauchy's theorem 
implies that the integral in Eq. ( 4.48j ) may be reduced to a line integral over the real 
axis: 



gen 



d(j,M(fi)V 



(4.49) 



4 We impose this restriction because it enables us to derive the 1- and 2-fermion states analytically. 
Solutions that depend on v can also be considered, although in such cases it is not possible to proceed 
analytically. 
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where M(/x) is an arbitrary function |39[ 



The Euclidean Hamilton-Jacobi equation (|4.44j ) becomes 

dldl , t . 

TsTr = 1 (4 ' 5 °) 

and admits the solutions 

I = _6s - _r. (4.51) 

where 6 = z/i. Eq. (|4.51|) is invariant under the duality transformation ( |2.12| ). More- 
over, we see the exact solution ( |4.47|) to the Wheeler-DeWitt equation ( |4.43| ) is also 
a WKB solution, \l/ = exp(±J), to the Euclidean Hamilton- Jacobi equation (|4.44|) . 

In performing the supersymmetric quantization of this cosmology, it is convenient 
to diagonalise the minisuperspace metric. The reason is that the Grassmannian vari- 
ables should satisfy the anticommuting relations [ip a , ip b ]+ = G ab . A non-diagonal 
minisuperspace metric would mean that fermionic states could not be clearly sepa- 
rated after the wavefunction has been annihilated by the supercharges. We therefore 
introduce the pair of variables 

T=i( a + T), X= l -{s-r) (4.52) 

and this implies that the minisuperspace metric, defined in Eq. ( [2.22 ), has the non- 
trivial components G 00 = -G 11 = -(A 2 - B 2 ){T 2 - X 2 )/2 and G 22 = 1. 

The supercharge ( p. 28 ) and its Hermitian conjugate (|2.29|) are then given by 

and 

q = _ijo^ _ mi JL _ ie 2 ^- - ie ^- - ie 1 ^-, (4.54) 

^ OT OX dv dT dX' v ; 

respectively, where 

6 + s) T -(»-*)*■ (4 - 55) 

For the supersymmetric wavefunction, we consider the ansatz 

= a+ + (3 b 6 b + -e abc -f c a b + a^6 1 6 2 , (4.56) 
2 

where {a±, (3b, are bosonic functions of the minisuperspace variables. The anni- 
hilation of the wavefunction ( (4.56j ) by the supercharges (|4.53 ) and ( [4.54| ) yields the 
set of coupled, first-order partial differential equations 

G °°w + Gll §§ + G22 fr + G °° { b + l) (3 °- ° u {~ b + l) (3l = (4 - 57) 
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G^-G^ + G^ (b + l T )^ 



dT dv \ b 



~ «F + 8X~ V + 6 J 7 G " l" 6 + 6/ 1 



da_ 
dv 







da- (. V. 

+ [b + ~ a_ = 



<9T V 6 

+ (6 - T a_ = 



<9X V b 
9a+ ^ + i , )o + = 



dT V 6 
<9a + / 1 



<9X V 6 



a + = 
<9a;+ 





<9t> 



4.59) 
4.60) 
4.61) 
4.62) 
4.63) 
4.64) 
4.65) 
4.66) 
4.67) 
4.68) 
4.69) 
,4.70) 

dT dX dv V b) 1 V b) 1 v ; 

The wavefunctions for the empty and filled fermion sectors are readily deduced: 

a± = e Tl , (4.71) 

where I is given by Eq. ( f4.55| ). To solve the remaining equations, we assume that the 
amplitudes {/3&, j c } are independent of the variable, v. In this case, Eqs. ( |4.68|) and 
( 14.691 ) yield the general solution, /3 2 = exp(— 7), modulo a constant of proportionality, 
where I is given by the Euclidean action ( [4.511) . Likewise, Eqs. fl4.58|) and (|4.59|) imply 
that 7 2 = exp(J). 

The wavefunctions for the one-fermion sector are completely determined by solv- 
ing Eqs. ( J4.57| ) an d ( 14.671 ). To proceed, it is convenient to transform back to the null 
coordinate pair (s, r) defined in Eq. (|4.45 ). In terms of these variables, we find that 

d W° + M +b{ p -{ k) = Q . (4.73) 



dT dX V 6/ V b 



/3 2 = 



dT dv V 6 

d/3 a / i 



«9X V b 
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Defining Y = f3 ~ A and Z = 6 _1 (/? + Pi) implies that Eqs. ( $~7%) and (|473D may 
be expressed in the more compact form: 



ds ~ J dr 



Y. 



(4.74) 



Differentiating the first constraint in Eq. ( 4.74j) with respect to r and substituting 
in the second condition implies that both /?o,i satisfy the unit-mass Klein-Gordon 
equation, i.e., the bosonic Wheeler-DeWitt equation ( [4.46| ). Thus, although these 
amplitudes satisfy the same equation as the wavefunction that arises in the standard 
quantum cosmological approach, the supersymmetry imposes strong constraints, as 
summarized in Eq. Q4.74j ), on the functional form of the solutions that can arise. One 
class of allowed solution is given by Y = exp(-I) and Z = —bY, where / is given by 
Eq. (pip. 

It now only remains to solve Eqs. (|4.60|) and (|4.70| ) in order to determine the two- 
fermion sector of the supersymmetric wavefunction. By combining and subtracting 
these two equations, we find that 



g(7° + 7 1 ) 
ds 

dr 



b ( 7 ° + 7 1 ) = 







0. 



(4.75) 
(4.76) 
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) implies that Eqs. ( 4.75|) and (|4.76 ) are 
dW 



R. 



(4.77) 



Defining R = 7 + 7 1 and W = b - 
equivalent to: 

dR 

— = W, 
ds dr 

Thus, the amplitudes 7 ' 1 also satisfy the unit-mass Klein-Gordon equation. We find 
that one class of solution consistent with Eq. ( [4.77D is given by R = exp(7) and 
W = —bR. To summarize, therefore, the supersymmetric wavefunction that we have 
found for the vacuum Brans-Dicke, Kantowski-Sachs cosmology is given by 



(4.78) 



where /3o,i and 7 ' 1 satisfy the unit-mass, Klein-Gordon equation (|4.46|) subject to 
the integrability conditions (|4 . 74|) and (|4.77|). 



5 Discussion 

In this paper, we have considered an N = 2 supersymmetric quantization of the 
LRS Bianchi type I and Kantowski-Sachs, Brans-Dicke cosmologies. In the former 
case, we found that a supersymmetric quantization is possible if a negative cosmo- 
logical constant is introduced into the gravitational sector of the theory. For the 
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Kantowski-Sachs universe, the existence of such a term is not necessary, because 
this model has positive spatial curvature. In both models, supersymmetric quantum 
states were found for a given solution to the Euclidean Hamilton- Jacobi equation. 
Furthermore, these wavef unctions respect a global scale factor duality symmetry of 
the respective classical Hamiltonians. 

Having found particular solutions to the supersymmetric quantum constraints, 
the immediate question that arises is the nature of the boundary conditions that such 
solutions satisfy. In general, the supersymmetric Hamiltonian has a spin term with 
a coefficient determined by the solution to the Euclidean Hamilton-Jacobi equation. 
This term implies that it is very difficult to complete a supersymmetric extension of 
the system with complex or imaginary solutions |17|]. Thus, the boundary conditions 
that are typically most relevant in this quantization scheme are those due to Hartle 



and Hawking [19] and to Hawking and Page [21, 22 



In particular, it is natural to consider whether the Kantowski-Sachs wavefunction 
derived above satisfies the Hawking-Page boundary conditions relevant to a wormhole 
configuration ETI, B^]. Classically, a wormhole represents an instanton solution of 



the Euclidean field equations |4(J |2T|. At the quantum level, such a state may be 
interpreted as a solution to the Wheeler-DeWitt equation. The appropriate boundary 
conditions that must be satisfied are that the wavefunction should be regular, in the 
sense that it does not oscillate an infinite number of times, when the three-metric 
degenerates and that it should be exponentially damped when the three-geometry 



tends to infinity plL 22 



The anisotropic geometry, S 1 x S 2 , of the Kantowski-Sachs model implies that 
there are different types of possible wormholes |4l], |42|. These have been studied by 
Campbell and Garay within the context of Einstein gravity minimally coupled to a 
massless scalar field |42| . The geometry of the spacetime asymptotes to R 3 x S 1 if the 
radius of the circle, Si, tends to a constant as the radius of the two-sphere diverges. 
Alternatively, if the volume of the two-sphere tends to a constant as di — > oo, the 
geometry is R 2 x S 2 . The wavefunction representing the ground state of each of 
these wormholes is the path integral over all metrics that asymptotically have these 
geometries and over all matter configurations that vanish at infinity. For the R 3 x S 1 
wormhole, the wavefunction is given by \& oc exp(— Adid 2 ) in the asymptotic limit. 
The corresponding limit for the R 2 x S 2 wormhole is \I/ oc exp(— a 2 ). 

After transforming back to the original variables of Section II, we find that the 
bosonic component of the supersymmetric Kantowski-Sachs wavefunction ( 4.78|) does 
not asymptote to either of these forms. Its interpretation as a quantum wormhole is 
therefore not clear. However, a further solution to the Euclidean Hamilton-Jacobi 
equation ( |4.44|) that respects the scale factor duality ( |2.12|) of the classical action is 
given by 

1 e^+ Bu )l\ (5.79) 



.A 2 - B 2 , 

Consequently, a supersymmetric quantization may be performed with this solution. 
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Due to the non-trivial functional form of Eq. ( |5.79|) , however, it has not been possible 
to find analytical solutions for the intermediate fermionic sectors. On the other hand, 
the empty fermion sector is given by \1/ oc e _/ and it is of interest to compare this 
wavefunction with the above ground state wormhole wavef unctions. For example, in 
the superstring inspired model, where u = —1, we find that / = 4aia2e~*. Performing 
a conformal transformation on the four-metric, g^ = Q 2 g fiU , where G 2 = e~*, implies 
that the dilaton field is minimally coupled to gravity in the 'Einstein-frame', g^. In 
terms of variables defined in this frame, the wavefunction is given by \l/ oc exp(— Aa^) 
and this is precisely the wavefunction for the R 3 x S 1 quantum wormhole that arises 
in the standard Wheeler-DeWitt quantization. 

This is important because the ground state of the R 3 x S 1 quantum wormhole has 
been selected by the supersymmetric quantization procedure. We emphasize that the 
bosonic component of the supersymmetric wavefunction is unique once a solution to 
the Euclidean Hamilton- Jacobi equation has been specified. In this sense, therefore, 
any ambiguities that arise in the operator ordering are eliminated. 

Moreover, the interior of a Euclidean Schwarzschild black hole has the form of a 
Kantowski-Sachs metric [^TJ and it is possible, therefore, that supersymmetric quan- 
tum cosmology may relate a black hole interior to a quantum wormhole. It would 
be interesting to consider this possibility further. For example, such a relationship 
would have implications for the graceful exit problem of the pre-big bang inflationary 
scenario P2[ . This problem arises because the classical, dilaton-driven inflationary 
solution becomes singular in a finite proper time. At present, no generally accepted 
mechanism has been proposed to avoid such a singularity and ensure a smooth tran- 
sition to the standard, post-big bang expansion. However, an epoch of pre-big bang 
inflation may be formally interpreted in the Einstein-frame in terms of gravitational 
collapse [[13] . If the final state of such a collapse were a non-singular supersymmetric 



wormhole configuration, such a problem could in principle be avoided. It is intrigu- 
ing that whereas the pre- and post-big bang branches are related to one another by 
the scale factor duality of the classical action, the empty fermion component of the 
wavefunction is invariant under such a duality transformation. 

In principle, the supersymmetric quantization of other homogeneous, scalar-tensor 
cosmologies can also be considered following the method outlined in this paper. The 
Bianchi type II, VIo and VIIo cosmologies also exhibit global symmetries at the clas- 
sical level and, in particular, the Wheeler-DeWitt equation for the Bianchi type 
II model reduces to Eq. ([4.46 ) after appropriate field redefinitions |^!J. Thus, a 
similar analysis to that presented in Section IV may also be performed for this 
model. Similarly, the effective potential arising in the Wheeler-DeWitt equation 
of the LRS type III model has an opposite sign to that given in Eq. ( f4.43|) . However, 
the Wheeler-DeWitt equation can be transformed into the unit-mass Klein-Gordon 
equation ( 4.46| ) after a suitable choice of null variables. 

Finally, there remains the open question of the possible relationship between the 
different approaches to supersymmetric quantum cosmology. As we discussed in the 
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introduction, a supersymmetric minisuperspace may be obtained directly from a full 
four-dimensional N — 1 supergravity action with the assistance of a suitable di- 
mensional reduction for both the bosonic and fermionic variables (see, e.g., Refs. 
H [U], [16]]). Alternatively, a bosonic minisuperspace may be extracted from a (1 + 0)- 
dimensional lagrangian and a supersymmetric extension established along the lines of 
0,0 



Refs. 24 



or 



27, 



Determining the fundamental similarities and differences 
between these two methods is a challenging problem. The one attempt to investigate 
this was made in Ref. 



but unfortunately it was based on an incomplete ansatz 
for the supersymmetric Bianchi type IX model. This particular problem of the ansatz 
was eventually corrected 0, [K], [L6| but no further studies have been made. Such a 
complex investigation is beyond the objectives and scope of this paper, but it would 
be interesting to consider this topic further. 



Acknowledgements 

This work was supported by the qrants ESO/INF/1260/98, ESO/PRO/1258/98, CRUP-BC 
N^B/73/99 and CERN/P/FIS/15190/1999. JEL is supported by the Royal Society. PVM 
is thankful to QMW for hospitality where some of this work was completed. The authors 
are grateful to M. Cavaglia and N. Kaloper for helpful discussions. 



References 

[1] Halliwell J 1990 Proceedings of the Jerusalem Winter School: Quantum Cosmology and 
Baby Universes eds. S Coleman, J Hartle, T Piran and S Weinberg (Singapore: World 
Scientific) 

[2] DeWitt B S 1967 Phys. Rev. 160 1113 

Wheeler 1968 Battelle Rencontres (New York: Benjamin) 

[3] Kuchaf K V and Ryan M P 1989 Phys. Rev. D 40 3982 

[4] Kuchaf K V and Ryan M P, 1986 in Gravitational Collapse and Relativity eds. H Sato 
and T Nakamura (Singapore: World Scientific) 

[5] Kuchaf K V, to appear in MG9, Rome 2000 Proceedings 

[6] Sinha S and Hu B L 1991 Phys. Rev. D 44 1028 

Hu B L, Paz J P and Sinha S 1993 Directions in General Relativity eds. B L Hu, M P 
Ryan and C V Vishveshwara (Cambridge: Cambridge University Press) 

[7] Henneaux M 1980 Gen. Rel. Grav. 12 137 

[8] Halliwell J J 1993 Proceedings of the Thirteenth International Conference on General 
Relativity eds. R J Gleiser, C N Kozameh and O M Moreschi (Bristol: Institute of 
Physics) 

[9] Moniz P V 1996 Int. J. Mod. Phys. A 11 4321 



17 



[10] D'Eath P D 1996 Supersymmetric Quantum Cosmology (Cambridge: Cambridge Uni- 
versity Press) 

[11] van Nieuwenhuizen P 1981 Phys. Rep. 68 189 

Wess J and Bagger J 1992 Supersymmetry and Supergravity (Princeton: Princeton 
University Press) 

[12] Green M B, Schwarz J H, and Witten E 1987 Superstring Theory, Vols. I and II 
(Cambridge: Cambridge University Press) 

Polchinski J 1998 String Theory, Vols. I and II (Cambridge: Cambridge University 
Press) 

[13] Witten E 1995 Nucl. Phys. B 443 85 
[14] D'Eath P D 1984 Phys. Rev. D 29 2199 

[15] Macias A, Obregon O and Ryan M 1987 Class. Quantum Grav. 4 1477 
Macias A, Obregon O and Socorro J 1993 Int. J. Mod. Phys. A 8 4291 

[16] Cheng A D Y, D'Eath P D and Moniz P V 1994 Phys. Rev. D 49 5246 

Cheng A D Y, D'Eath P D and Moniz P V 1995 Class. Quantum Grav. 12 1343 

Csordas A and Graham R 1995 Phys. Rev. Lett. 74 4129 

Cheng A D Y and Moniz P V 1995 Int. J. Mod. Phys. D 4 189 

Cheng A D Y and Moniz P V 1996 Mod. Phys. Lett. A 11 227 

Moniz P V 1996 Int. J. Mod. Phys. A 11 1763 

Moniz P V 1998 Phys. Rev. D 57 R7071 

[17] Graham R 1991 Phys. Rev. Lett. 67 1381 

Bene J and Graham R 1994 Phys. Rev. D 49 799 

[18] Lidsey J E 1995 Phys. Rev. D 51 6829 
Lidsey J E 1995 Phys. Rev. D 52 R5407 

[19] Hartle J and Hawking S W 1983 Phys. Rev. D 28 2960 

[20] Vilenkin A 1982 Phys. Lett. 117B 25 

[21] Hawking S W 1988 Phys. Rev. D 37 904 

[22] Hawking S W 1990 Mod. Phys. Lett. A 5 453 

Hawking S W and Page D N 1990 Phys. Rev. D 42 2655 

[23] Teitelboim C 1977 Phys. Rev. Lett. 38 1106 
Teitelboim C 1977 Phys. Lett. 69B 240 
Pilati M 1978 Nucl. Phys. B 132 138 

[24] Witten E 1981 Nucl. Phys. B 188 513 

[25] Claudson M and Halpern M B 1985 Nucl. Phys. B 250 689 
Graham R and Roeckaerts D 1986 Phys. Rev. D 34 2312 



18 



[26] Davies A C, MacFarlene A J, Popat P C, and van Holten J W 1984 J. Phys. A 17 
2945 



[27] Gozzi E 1983 Phys. Lett. 129B 432 

[28] Andrianov A A, Borisov N V, and Ioffe M V 1984 Phys. Lett. A 105 19 

Castanos O, D'Olivo J C, Hojman S, and Urrutia L F 1986 Phys. Lett. 174B 307 

Gozzi E 1983 Phys. Rev. D 28 1922 

Gamboa J and Zanelli J 1985 Phys. Lett. 165B 91 

Gamboa J and Zanelli J 1988 Ann. Phys. 188 239 

[29] Ryan M P and Shepley L C 1975 Homogeneous Relativistic Cosmologies (Princeton: 
Princeton University Press) 

MacCallum M A H 1979 General Relativity: An Einstein Centenary Survey eds. S W 
Hawking and W Israel (Cambridge: Cambridge University Press) 

[30] Brans C and Dicke R H 1961 Phys. Rev. 124 925 

[31] Veneziano G 1991 Phys. Lett. 265B 287 

[32] Gasperini M and Veneziano G 1993 Astropart. Phys. 1 317 



[33] Lidsey J E, Wands D and Copeland E J 2000 [hep-th/9909061 
Gasperini M and Veneziano G 2000 in preparation 

[34] Giveon A, Porrati M and Rabinovici E 1994 Phys. Rep. 244 77 

[35] Lidsey J E 1994 Phys. Rev. D 49 R599 

Bento M and Bertolami O 1994 Class. Quantum Grav. 11 1211 

Gasperini M, Maharana J, and Veneziano G 1996 Nucl. Phys. B 472 344 

Gasperini M and Veneziano G 1996 Gen. Rel. Grav. 28 130 

Kehagias A and Lukas A 1996 Nucl. Phys. B 477 549 

Dabrowski M and Kiefer C 1997 Phys. Lett. 397B 185 

Cavaglia M and de Alfaro V 1997 Gen. Rel. Grav. 29 773 

Lidsey J E 1997 Phys. Rev. D 55 3303 

Lukas A and Poppe R 1997 Mod. Phys. Lett. A 12 597 

Cavaglia M and Ungarelli C 1999 Class. Quantum Grav. 16 149 

[36] Lidsey J E 1996 Class. Quantum Grav. 13 2449 

Clancy D, Lidsey J E, and Tavakol R 1998 Class. Quantum Grav. 15 257 

[37] Di Pietra E and Demaret J 1999 Int. J. Mod. Phys. D 8 349 

[38] Freund P G O 1982 Nucl. Phys. B 209 146 

[39] Page D N 1991 J. Math. Phys. 32 3427 

[40] Giddings S B and Strominger A 1990 Nucl. Phys. B 306 890 

Halliwell J J and Laflamme R 1989 Class. Quantum Grav. 6 1839 



19 



[41] Keay B J and Laflamme R 1989 Phys. Rev. D 40 2118 

Cavaglia M, de Alfaro V, and de Felice F 1994 Phys. Rev. D 49 6493 
Cavaglia M 1994 Mod. Phys. Lett. A 9 1897 

[42] Campbell L M and Garay L J 1991 Phys. Lett. 254B 49 

[43] Buonanno A, Damour T, and Veneziano G 1999 Nucl. Phys. B 543 1275 
Feinstein A and Vazquez-Mozo M A 1998 Phys. Lett. 441B 40 
Feinstein A, Kunze K E, and Vazquez-Mozo M A 2000 |hep-th/0004094 

[44] Graham R 1993 Phys. Rev. D 48 1602 

[45] Lidsey J E 1995 Phys. Lett. 352B 207 



20 



